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We introduce for SU(2) quantum spin systems the Valence Bond Entanglement Entropy as a 
counting of valence bond spin singlets shared by two subsystems. For a large class of antiferromag- 
netic systems, it can be calculated in all dimensions with Quantum Monte Carlo simulations in the 
valence bond basis. We show numerically that this quantity displays all features of the von Neumann 
entanglement entropy for several one-dimensional systems. For two-dimensional Heisenberg models, 
we find a strict area law for a Valence Bond Solid state and multiplicative logarithmic corrections 
for the Neel phase. 
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Entanglement is probably the most prominent feature 
that allows to distinguish quantum from classical sys- 
tems. While being an inherent ingredient of quantum 
information processing, the study of entanglement prop- 
erties of quantum many-body states has recently shed 
new light on our understanding of condensed matter sys- 
tems at zero temperature Over the past few years, 
various measures of entanglement have been used to in- 
vestigate quantum phase transitions and states of mat- 
ter P, [3, H 3 • Amongst them, the von Neumann entan- 
glement entropy (EE) quantifies the bipartite entangle- 
ment between two parts of a quantum system [Bj]. The 
EE of a quantum state \^>) between a part O and the rest 
of a system is 

Sn = - Trpn lnpn 

where pn = Tiy 2 |^ r )( 1 J r | is the reduced density matrix of 
£1 obtained by tracing out the rest of the system f2. An 
important property of EE is that SV2 — and it ap- 
pears naturally that Sq is only related to the common 
property of O and Q, their boundary: general arguments 
indeed indicate that Sq typically scales with the size of 
the boundary (area law) (gj. For critical systems how- 
ever, logarithmic corrections can be present. If the crit- 
ical system is conformal invariant, the amplitude of the 
logarithmic corrections is related to the central charge 
of the corresponding conformal field theory (CFT). This 
has been shown in one dimension (Id) |3|, [Z| and more 
recently for some CFT in two dimensions (2d) [H, in 
which case the coefficient also depends on the geometry 
of Q. For systems possessing topological order, other sub- 
leading corrections also depend on the topology of Q . 
The scaling of the EE with the size of therefore con- 
tains precious informations on the state of the physical 
system and can be used e.g. to detect criticality. This 
has been successfully demonstrated in quantum spin sys- 
tems 0, d, 0, EH El > which are natural candidates for 
such studies as they are both of theoretical and experi- 
mental relevance in condensed matter. 

While many interesting properties of EE have been 




FIG. 1: (Color online) A bipartite valence bond state. The 
valence bond entanglement entropy of the area f2 (shaded 
square) is equal to the number of singlets shared between Q 
and the rest of the system Q (here 8), times ln(2). 



derived exactly for intcgrable models or systems with ex- 
actly known ground-states (GS), the calculation of EE for 
general interacting quantum systems is an exacting task. 
On the practical side, even numerical simulations are dif- 
ficult. Exact diagonalization is limited to small systems 
and cannot precisely verify scaling properties and EE is 
not accessible to Quantum Monte Carlo (QMC) meth- 
ods. In Id, EE can be calculated within the DMRG 
method [l2| but this technique is no longer available in 
higher dimensions. In this Letter, we introduce another 
measure of entanglement, called Valence Bond Entan- 
glement Entropy (denoted hereafter S ), which is de- 
fined for quantum spin systems with SU(2) symmetry. 
It can be measured efficiently in all dimensions for non- 
frustrated antiferromagnetic (AF) spin systems via QMC 
simulations in the Valence Bond (VB) basis [l3| . We 
show that S VB , albeit different, captures all features of 
the von Neumann EE. In particular, it obeys the same 
scaling properties: for Id quantum spin systems, we re- 
cover CFT predictions. We also investigate the proper- 
ties of S VB for a 2d Heisenberg model possessing Neel and 
Valence Bond Solid phases. Our numerical results show 
that S VB obeys an area law for both phases, with the 
presence of logarithmic corrections for the Neel phase. 
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We focus on the AF spin-1/2 Heisenbcrg Hamiltonian 



H 



7 , Jij &i 



which conserves the total spin St of the system. AF 
interactions Jij > favor a singlet St = ground-state. 

Definition — It is well-known that any singlet state 
can be expressed in the VB basis, where spins couple 
pairwise in singlets (| fj.) — | J.T))/v2. The VB basis is 
overcomplete (there are more VB coverings than the to- 
tal number of singlets). Another basis is the bipartite VB 
basis [l4[ , where the system is decomposed into two sets 
(for instance different sublattices for a bipartite lattice) 
such that two spins forming a singlet necessarily belong 
to different sets. This basis is smaller, but still overcom- 
plete. In a given bipartition of sites and bipartite VB 
state \^f) such as the ones in Fig. [TJ consider a subsys- 
tem ft (shaded area in the figure) . We define the Valence 
Bond Entanglement Entropy of this state as: 

SX B (l*»=m(2). „<-(|*» 

where n^d^)) is the number of singlets that cross the 
boundary of Q (i.e. singlets with one and only one con- 
stituting spin within ft). The constant ln(2) is used to 
match the EE for a unique singlet. In general, the GS of 
the Hamiltonian H is not a single VB state. For a lin- 
ear combination |$) = J^i^l^i) with bipartite VB 
states, we define S^ B (|$)) = Ei ^cPCI**))/ Ei An 
estimate of S VB for the GS of H can then be obtained 
by using the Projector QMC method recently proposed 
by Sandvik (l3j , which precisely works in a bipartite VB 
basis (as soon as H is not frustrated). In this method, 
the GS is sampled by applying H n stochastically (with 
n sufficiently large) to an initial VB state |\&o) m order 
to obtain a projected state \*$> n ). In this process, \^f n ) 
appears a number of times proportional to its coefficient 
in the GS wave function. The GS estimate of S VB is 
obtained by Monte Carlo-averaging S' VB (|4 , „)) over the 
projected states obtained in all Monte Carlo steps. 

Bipartite basis and overcompleteness At first 

glance, S VB seems ill-defined as the bipartite VB basis 
is overcomplete. |$) can indeed be rewritten as a lin- 
ear combination of other states for which S VB could be 
different. However, S VB turns out to be conserved in 
any linear combination between the bipartite VB states. 
Via an algebraic representation of the singlet wave func- 
tion, it can be indeed shown that for any linear relation 
EjCil^i) = between the bipartite VB states we 
have E l c *< S fi B (|1'*)) = for all ft Another issue 

is the choice of the bipartite basis: in the general case, 
S VB does depend on the precise choice of bipartition. 
However, if |$) satisfies a Marshall sign criterion a 
genuine bipartition where all a% > exists and should be 
taken. Not all singlet states satisfy this, but this is the 
case for e.g. the GS of H on bipartite lattices In 



practice, the QMC method is efficient only for bipartite 
H (when all a, > 0) and the Marshall condition is not a 
restriction. This implies that S VB is a well-defined quan- 
tity for any singlet wave function satisfying the Marshall 
criterion and that it can easily be measured (e.g. graph- 
ically as in Fig. [T]) for any bipartite VB state. 

Entanglement properties — S VB is in general differ- 
ent from EE (this can be seen by computing exactly both 
quantities on a small system). However, S VB share many 
properties with EE: (i) S is a sub-additive quantity: 
S%% 2 < So? + Sff, (ii) S^ B = 5V B as any singlet 
crossing the boundary belongs to both subsystems, (iii) 
5^ B = Sq if ft contains a single site, (iv) Sj( B = Sq for 
all ft when the GS is a single VB state: this is the case in 
the thermodynamic limit for the Majumdar- Ghosh spin 
chain 17 1 or for the Random Singlet phase [l8| of disor- 
dered spin chains where one VB state essentially domi- 
nates |4|. Finally, we note that S VB offers a simple ge- 
ometrical interpretation of entanglement properties of a 
quantum spin state: ft is entangled with the rest of the 
system ft iff there are singlets in between these two parts. 

In the rest of this paper, we will show that S VB cap- 
tures both qualitative and quantitative features of the EE. 
S VB can be calculated for all systems that can be simu- 
lated with the VB QMC method, in particular all spin- 
1/2 AF Heisenberg-like models on bipartite lattices in all 
dimensions, including models with multiple-spin interac- 
tions [Tij . We first start with Id models, where all the 
characteristics of the EE are recovered with high preci- 
sion. We then present results on entanglement properties 
of 2d AF Heisenbcrg models. 

Id systems — We first consider Id uniform AF Heisen- 
berg chains where the GS is known to be critical, with 
algebraic decay of spin correlations. For systems of fi- 
nite size L with periodic boundary conditions (PBC), 
CFT predicts that the EE of a block of spins of size 
x scales (for large a;) as S(x) = c/31n(a;') + Sq, where 
x' = L/TTsm(nx/L) is the conformal distance, the cen- 
tral charge c = 1 and So a non-universal constant 0]. 
Fig. displays S VB (a/) for a Id uniform chain of size 
L = 128 and the resulting curve is well fitted by the 
form S VB (x') = c off /31n(a;') + S with c off = 1, in full 
agreement with CFT predictions for EE. 

Now we turn to uniform AF Heisenbcrg chains with 
open boundary conditions (OBC). For a segment of size 
x starting at the open boundary, the CFT prediction for 
EE is S(x) = c/61n(x') + Si with c = 1 and Si another 
constant related to So Our results (see Fig. (2}d) for 
S VB (x) for an L = 128 open chain also clearly follow this 
scaling form. We also recover the alternating term found 
for the EE of open Heisenbcrg chains [ll| , as can be seen 
from the two distinct odd-even sets of points in Fig. [2b- 

We next consider an alternating Heisenberg spin chain, 
where the coupling is 1 + 5 (respectively 1 — 6) for all even 
(resp. odd) bonds. This system is non-critical with a fi- 
nite correlation length for all 6 ^ and the EE S(x) of 
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FIG. 2: (Color online) 5 VB versus distance for Id AF Heisen- 
berg quantum spin systems: (a) chain with PBC, (b) chain 
with OBC, (c) alternating chain (5 = 1/3) with PBC, (d) 
random AF chain with PBC. In (a), (b) and (d), scale is log- 
linear. For (a) and (b), x 1 = L/n sin(irx/L) is the conformal 
distance. Solid lines are scaling forms for EE predicted by 
CFT [(a) and (b)] and RSRG [(d)]. Simulation parameters 
are : (a) n = 70L, (b) n = 60L, (c) n = 30L, (d) n = 10L 
and 500 disorder samples. 



a segment of size x saturates to a finite value for large 
enough x @, We calculate Syb{x) for 5 = 1/3 and a 
system of size L = 64 (see Fig. [2]:). There are two types 
of inequivalent blocks x depending on whether the block 
starts with a strong bond 1+5 or a weak one 1—5. We de- 
note the corresponding entropies by Sj B (x) and S^ B (x) 
and by symmetry Sj B (x) = S^ B (x) for all odd x. Both 
entropies show an alternating pattern, a consequence of 
the explicit dimerization. This is in contrast with the 
average entropy S VB (x) = l/2[Sj B {x) + SZ B {x)} which 
is a smooth function of x. All entropies clearly display a 
saturation at large x. These features also appear in the 
study of EE for the same model. 

We conclude this section by showing results on the 
disordered AF Heisenberg chain. Here all couplings are 
random variables uniformly distributed i n [ 0, 1], This 
system is in the Random Singlet phase [18[, a gapless 
critical phase where spins are coupled pairwise in sin- 
glets at all length scales. As understood from a Real 
Space Renormalization Group (RSRG) procedure [if] ], 
the physics is dominated by a single VB covering in the 
thermodynamic limit fl8j . Calculating the disorder aver- 
age EE S(x) wit h RSR G, Refael and Moore 0] obtained 
the scaling form S(x) = 7/31n(.-r) + 5 2 with 7 = ln(2), a 
result confirmed by numerical simulations of random spin 
chains [2l| . We simulated disordered chains of sizes up to 
L = 256 starting from an initial VB state l^o) obtained 
applying the RSRG decimation scheme, which ensured a 



faster convergence with the projection index n. Fig. [2ji 
presents our numerical results for S VB (x), and the data 
(for 10 < x < L/2) are well-fitted by the previous func- 
tional form with a prcfactor 7 c ft — 0.6(1), in agreement 
with the RSRG prediction for the average EE. Actually, 
as the average in the thermodynamic limit is dominated 
by a single VB state, we expect 5 VB and S to coincide. 
This is however difficult to check numerically because of 
the finite L and the small number of random samples used 
in the simulations. This also accounts for the small dis- 
crepancy in our numerical estimate of 7 c ff with respect to 
the RSRG prediction (note that a fit forcing 7 e g = ln(2) 
is also of good quality as seen in Fig. [2Ji). 

2d systems — The previous section showed the precise 
agreement between the VB entanglement entropy S VB 
and the true EE S in both quantitative and qualitative 
respects for Id systems. This gives us confidence that 
5 VB is a good observable to quantify entanglement also 
in larger dimensions d. This is an important point as 
the VB QMC simulations (and therefore the calculation 
of 5 VB ) can easily be extended to higher d. Note at this 
stage that only few results are available on entanglement 
properties in high d: exact results have been derived for 
free fermions 22J, 23[, bosons @,[23, 24 1 and some exactly 
solvable spin models [25j], but EE has never been calcu- 
lated numerically for large d > 1 systems to our best 
knowledge [26j . At present time, the connection between 
the area law and the precise nature of the ground-state 
(range of correlations, presence of a gap) is not fully un- 
derstood for d > 1 (see a discussion in Ref. 1241a). 

We now investigate 2d Heisenberg systems and calcu- 
late S VB for a square subsystem of linear size x (such 
as the one in Fig. [IJ, for L X L samples with PBC up 
to L = 64. In particular, we consider a model of cou- 
pled dimers in two dimensions depicted in the inset of 
Fig. [3j where dimers (with intra-dimer exchange J = 1) 
are coupled with an inter-dimer exchange A. The physics 
of this model is well-understood [27j : at low A, the sys- 
tem is in a dimerized Valence Bond Solid (VBS) phase 
whereas the case A = 1 corresponds to the Neel phase 
of the isotropic 2d Heisenberg model. The system un- 
dergoes a quantum phase transition between these two 
states at A c = 0.52337(3) [22j|. Fig. [3K displays the aver- 
age Valence Bond entanglement entropy divided by the 
linear size S VB (x)/x for the point A = 0.2 located in the 
VBS phase. All the curves for different L saturate for 
large enough x, indicating that the Valence Bond entan- 
glement entropy obeys an area law S YB (x) oc x, as ex- 
pected for a gapped phase. For the isotropic case A = 1 
located in the Neel phase, S VB (x)/x does not saturate, 
and the curve is actually well-fitted by a log form, as 
exemplified by the log-linear scale of Fig. (3}d. Our nu- 
merical results therefore indicate that the Valence Bond 
entanglement entropy scales as 5 VB (cc) oc x\n(x) in the 
Neel phase. Multiplicative logarithmic corrections have 
also been found for the EE of fermions 0, [23|. They 
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FIG. 3: (Color online) S VB /x versus distance x for 2d quan- 
tum spin systems: (a) 2d AF dimerized Heisenberg model 
(A = 0.2), (b) 2d AF Heisenberg model (A = 1) (log-linear 
scale). Solid line is a fit to a logarithmic divergence form. 
Inset: Illustration of the coupled-dimer model. Simulation 
parameters are n = 30L 2 for L < 32, n = 20L 2 for L > 32. 



are tentatively attributed here to the Goldstonc modes 
present in the Neel phase, which also make for the alge- 
braic transverse spin correlations. Scaling of S YB at the 
critical point A c is currently under investigation. 

To summarize, we have proposed the Valence Bond en- 
tanglement entropy as a new measure of entanglement for 
SU(2) quantum spin systems. S VB offers a powerful tool 
to quantitatively investigate the deep connection between 
VB physics and the amount of entanglement present in 
spin systems, although its precise connection to EE needs 
further studies. Several questions then arise, pointing 
for instance towards d > 2 exotic phases (such as RVB 
spin liquids) and quantum critical points (as in Ref. [3a), 
which should hopefully exhibit striking valence bond en- 
tanglement features. Finally, in the absence of a Marshall 
sign rule (as for some frustrated systems) where S VB de- 
pends on the choice of a bipartite basis, the usefulness of 
S YB has to be clarified. 
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